Low-Energy Linear Structures in Dense Oxygen: Implications for the e-phase 
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Using density functional theory implemented within the generalized gradient approximation, a new 
non-magnetic insulating ground state of solid oxygen is proposed and found to be energetically fa- 
vored at pressures corresponding to the e-phase. The newly-predicted ground state is composed of 
linear herringbone-type chains of O2 molecules and has Cmcm symmetry (with an alternative mon- 
oclinic cell). Importantly, this phase supports IR-active zone-center phonons, and their computed 
frequencies are found to be in broad agreement with recent infrared absorption experiments. 



PACS: 61.50.Ah, 78.30.Am 

At low temperatures and ordinary pressures solid oxy- 
gen condenses into the only antiferromagnetic insulating 
phase known among the elemental solids. In its ground 
state, the well-known monoclinic a-phase Q, the molec- 
ular spins are arranged in nearly-closed packed planes (or 
layers) perpendicular to the molecular axes; the axes of 
molecules in successive planes are collinear |lj-|3|] (and ap- 
parently so in all measured phases) . As the temperature 
is increased, oxygen eventually undergoes a transition to 
its familiar room-temperature gaseous phase, but as the 
pressure is increased, its attributes depart radically from 
those of an ensemble of weakly-interacting molecules. In 
fact at pressures above 96 GPa (at nearly 3-fold compres- 
sion) , both diamond-anvil and shock j|] experiments 
have reported an insulator-metal transition and, at tem- 
peratures below 0.6 K, the metallic state (the £-phase) 
even exhibits superconductivity While some of the 
broader aspects of the electronic behavior are understood 
throughout this density range, the crystal structure and 
magnetic properties of the so-called e-phase, persisting 
as it does over the wide, intermediate pressure range of 
10-96 GPa, remain unknown. 

Over two decades ago, a dramatic color change was ob- 
served (from light-blue transparent to darkening red) in 
experiments above 10 GPa at room temperature ||. De- 
tailed optical measurements || shortly thereafter at room 
temperature revealed an abrupt onset of infrared (IR) ab- 
sorption in the fundamental molecular vibron, also at ap- 
proximately 10 GPa. Recent low-frequency optical stud- 
ies have uncovered another strong IR peak, ranging in 
frequency between 300-500 cm -1 from 10 to 70 GPa, but 
at very much lower frequencies than that of the vibron 
fTs] |l5| . The formation of pairs of O2 molecules with T)2h 
(rectangular) symmetry-an O4 unit-has been suggested 
to explain these new low- frequency modes |T^ , |l5|] . Inter- 
estingly enough the sudden increase in the intensity of 
infrared activity in two quite different frequency bands 
parallels the case of dense hydrogen, where a spontaneous 
polarization of the charge density along the molecular 



bonds has been suggested to account for the stability of 
low-symmetry structures above 150 GPa ]l~6[ ]. 

Although there exist several promising structures 
H |l4| , to date x-ray diffraction studies have been unable 
to completely reveal the exact positions of the molecules 
in the e-phase. Likewise a previous first-principles study 
recorded a premature magnetic collapse and concomitant 
metallization into to the £-phase near 12 GPa pTj , by- 
passing the e-phase. In this Letter we summarize results 
of first-principles calculations that predict a new low- 
enthalpy molecular arrangement in the pressure range of 
the e-phase. Beginning with a symmetrical low-density 
structure, we observe that it is unstable to the formation 
of extended herringbone-type chains of O2 molecules (in- 
stead of O4 units), strikingly similar to that suggested 
by Agnew et al. jnj. Non-magnetic and insulating, this 
newly-predicted phase is also consistent with recent in- 
frared measurements. 

We examine here the stability of the common symmet- 
ric phases of oxygen over the molar volume range 6.7- 
13 cm 3 /mole, the lower value corresponding to pressures 
exceeding 100 GPa. Density functional theory (DFT) is 
used within the local spin-density approximation (LSDA) 
]l8[ and with gradient corrections p9| . We utilize the 
projector augmented-wave (PAW) method [^(J as imple- 
mented within the Vienna ab initio Simulations Package 
(VASP) Q. Our oxygen PAW potential relegates the Is 
electrons to a frozen core but otherwise treats all other 
electrons explicitly as valence; a 60 Ry plane wave cut- 
off is used for all calculations. These methods provide 
a particularly accurate picture of the free molecule for 
which we correctly obtain the magnetic (5 = 1) ground 
state; our calculated bond length of 1.236 A, binding en- 
ergy of near 6.0 eV, and molecular vibron frequency of 
1550 cm -1 (within the harmonic approximation) are all 
slightly larger than experiment (1.207 A) but quite con- 
sistent with the known tendency of gradient corrections 
to overestimate the bond length and binding energy of 
p-bonded diatomic molecules p3]. 
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A primary physical issue centers on the notably high 
linear dipole polarizability of the oxygen molecule. The 
along-axis tensor component is 15.9q,q and the off-axis 
components are 8.2oq pH , where ao is the Bohr radius. 
These exceed the corresponding quantities in hydrogen 
by more than a factor of two, and the Hertzfeld crite- 
rion (for the onset of a polarization divergence) would 
require a compression of only about 1.8. There is a 
large anisotropic fluctuating-dipole (or van der Waals) 
attraction [Q, and an accurate, effective, and fully 
non-local representation of these correlated fluctuations 
within density functional theory has proven to be a chal- 
lenge |2^,^(| . Thus the local density approximation is not 
expected to be satisfactory until significant intermolecu- 
lar density has accumulated, and accordingly we focus 
our study on volumes less than 13 cm 3 /mole (8 GPa), 
above which we obtain adequate agreement (± 5 GPa) 
with the equation of state at 300 K [[|. 
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FIG. 1. Enthalpy (H — E + pV) vs. pressure p for selected 
structures of dense oxygen with respect to the Fmmm phase 
at T=0 K. In general agreement with a previous study Jl^ , we 
observe that the Fmmm phase is unstable to C2mm above 
17 GPa. Convergence with respect to k-points is achieved at 
1 meV/molecule using 10x10x12 Monkhorst-Pack k- meshes 
for 2-molecule primitive cells; more k-points are used for 
the metallic phase. Relaxations are performed until Hell- 
mann-Feynman forces are less than 10 -2 eV/A. The enthalpy 
H was determined by fitting the energy E to X^n=-2 a iV~ n ^ 3 ■ 

At room temperature and under a moderate com- 
pression of 8 GPa, the antiferromagnetic r5-phase (space 
group Fmmm) of solid oxygen has been particularly well- 
characterized [^HID . However recent x-ray data suggests 
a direct a-to-e transition at low temperatures Jl^| . Fur- 
ther, IR spectra @ and Raman j27) between 2-8 GPa 
and below 20 K appear to be inconsistent with Fmmm 
symmetry. Thus in order to examine the possibility of 
lower symmetry ground states that do possess IR activity 
at the observed frequencies and remain insulating, we re- 
leased the symmetry constraints of a 4-molecule Fmmm 
orthorhombic cell and completely relaxed the internal co- 
ordinates and lattice parameters, all with dense k-point 
sampling. This immediately resulted in a significant rear- 
rangement of the crystal and a new, stable orthorhombic 



structure, possessing a two-molecule primitive cell and 
Cmcm space-group symmetry. A plot of the enthalpies 
as a function of pressure for all structures considered ap- 
pears in Fig. 1. For comparison, we also evaluate the 
enthalpy of the metallic Clmm phase, a candidate 
phase proposed by Serra et al. [|l7| . Interestingly we find 
the Cmcm structure to be unstable to Clmm at a calcu- 
lated pressure of 47 GPa, resulting in an insulator-metal 
(IM) transition at a lower pressure than experimentally 
observed (96 GPa Q). Since the IM transition is associ- 
ated with band gap closure, this discrepancy may result 
from the usual underestimate of the electronic band gap 
by the gradient-corrected LDSA. 
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FIG. 2. Schematics of (010) ac-planes of (a) Fmmm, (b) 
C2mm, and (c) newly-predicted Cmcm phases. The molecules 
in white (oriented along [001]) are shifted by half a lattice vec- 
tor out of the plane. The orthorhombic Cmcm phase (having 
Wyckoff positions 8(g)) results from an intraplanar distortion 
in which the center molecule moves off-site. Chains run per- 
pendicular to the plane and are indicated by with bold gray 
lines. Differences in, e.g., a or c/a are ignored in these illus- 
trations. 

Atomic arrangements for the three phases considered 
here appear in Fig. 2. In the Cmcm structure the oxygen 
molecules order into symmetric herringbone-type chains 
along [010] (the b- axis) and perpendicular to the molec- 
ular bonds. They form through the shearing of adjacent 
(010) planes of the Fmmm structure, reducing the coor- 
dination of the molecules from four to two. An e-phase 
having monoclinic symmetry has been suggested j|,|],|l2| , 
and in this context we note that a monoclinic primi- 
tive lattice vector c' can be chosen for Cmcm (where 
|c'| = i Va 2 + c 2 ). (A similar relationship exists between 
the a (C2/m) and S (Fmmm) phases Jl2| .) Nevertheless 
we find that distortions of its orthorhombic lattice vectors 
a and c (resulting in P2\/c symmetry) do not lower the 
total energy f28f| . The Fmmm phase is related through 
continuous distortion to the nearly close-packed Clmm 
phase by a displacement of (001) planes in the [100] di- 
rection [^9|. Interestingly, the dimer length is found to 
decrease by about 1% with increasing pressure in this 
new phase, from 1.224 A at 8 GPa to 1.209 A at 54 GPa. 
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The distance between neighboring molecules along the 
chain is reduced from 2.081 A to 1.986 A over this range, 
a modest overbinding with respect to experiment, where 
these distances are thought to be in the range 2.2-2.5 A 
p3[ . Similarly, the a, 6, and c lattice parameters decline 
from 3.805 to 3.333 A, 2.990 to 2.718 A, and 7.034 to 
6.216 A, respectively, as pressure is increased from 8 to 
54 GPa. Notably, c diminishes slightly faster than either 
a or b with increasing pressure; experimentally, the &-axis 
is observed to decrease most rapidly |l4j. As the c-axis 
spacing is much larger than either a or b, remnant van der 
Waals interactions (see above) may still be particularly 
important between adjacent (001) planes, and therefore 
it is here that our treatment of exchange and correlation 
is likely to be most inadequate. 




FIG. 3. Charge density slice of Cmcm at a calculated 
pressure of 20 GPa (10.26 cm 3 /mole) plotted in an (001) (or 
ab) plane passing through the molecular centers (the unit cell 
is outlined in black). A natural logarithmic scale is used; the 
highest contour (blue) are 54 times larger than the lowest 
(red). Since the chains are equally spaced in the cell, the 
development of an alternating polarization perpendicular to 
the molecule is evident by inspection. 

Solid oxygen is robustly insulating in the Cmcm struc- 
ture, and the calculated band gap (a considerable un- 
derestimate of the true band gap) decreases from 0.95 
eV (direct gap) near 8 GPa to about 0.55 eV (indirect 
gap) at ~55 GPa. The optical threshold of the e-phase 
has been measured to be around 3 eV (blue) near 10 
GPa, declining to ~2.0 eV (orange-red) near 55 GPa @, 
for light polarized perpendicular to the molecular axis. 
These are larger than our computed direct gaps at these 
pressures (for the smallest, by a factor of three) but, 
as stated above, this is to be expected given our use of 
the LSDA; corrections through approximate inclusion of 
many-electron effects via, for example the GW approxi- 
mation |3(|, would be of considerable interest for accu- 
rately reproducing this gap and likewise describing the 
onset of metallization. 

The charge density is plotted in Fig. 3 in the a6-plane 
(perpendicular to [001]), and is seen to develop an asym- 
metry occurring perpendicular to the molecular bond but 
centered on each molecule, as is evident by inspection of 
Fig. 3, providing insight into the stability of this phase. 
The asymmetry may be viewed as a weakly antiferroelec- 



tric state, whose physical origin may in turn be under- 
stood through consideration of a mean-field argument for 
a dynamic lattice as given for hydrogen p6| . Thus self- 
consistent Lorenz fields compensate energetic costs as- 
sociated with the low-symmetry Cmcm structure in the 
intermediate pressure range below 50 GPa. Increasing 
intermolecular overlap at higher densities lowers kinetic 
and exchange energies, and the metallic C2mm phase is 
thus eventually favored. Examination of the spin density 
indicates that oxygen is non-magnetic within this phase, 
consistent with other studies p3| , pl| ; total energies and 
forces obtained from spin- and non-spin-polarized calcu- 
lations are essentially identical. 
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FIG. 4. Calculated zone-center phonons of Cmcm as a 
function of pressure. Modes having ~&2u and B3„ symmetry 
are IR-active (red-dashed); the remaining are Raman-active 
(black-solid). Lines guide the eye. 

Since there are two molecules in each primitive Cmcm 
cell, we expect nine optical phonons (plus three purely 
translational modes). At the zone-center, there are 
eight irreducible representations, as permitted by the 
D2/1 point group: four allow Raman activity (A s , Bi g , 
B 2g , and B lg ), three are IR-active (B lu , B 2u , and B 3u ), 
and a single remaining A u mode is silent. Gorelli et 
al. [ fb3|Jl5f| have proposed a structure made up of two 
molecules (actually O4 units) having T>2h symmetry. In 
agreement with Ref . |l5| , our predicted ground state pos- 
sesses T>2h point symmetry, but also includes additional 
lattice translations associated with the herringbone-type 
chain structure. Corresponding force constants are ob- 
tained through analysis |?2| of a series of frozen-phonon 
calculations; phonon frequencies are then computed by 
diagonalization of block-diagonal dynamical matrices. 

The pressure dependence of zone-center optical fre- 
quencies computed between 8 and 86 GPa appears in 
Fig. 4, and we obtain remarkably good agreement with 
the existing spectroscopic data, quite apart from the 
fact that C2mm is favored above 47 GPa in our cal- 
culations. Consistent with observations jl5| , we predict 
two IR-active modes over this range. The B3„ mode is 
an antiphase vibron; its computed frequencies are in the 
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range 1370-1480 cm -1 , ywhich, while slightly lower than 
found in experiment [fL5[ , increase with increasing pres- 
sure in qualitative accord with measurements and a de- 
creasing dimer length. The B2« mode is an antiphase 
libron-like mode with frequencies between 465-703 cm -1 
over the pressure range 8-86 GPa. As with the vibron, 
these frequencies are somewhat larger than those found 
experimentally, which increase from 300 to 600 cm -1 be- 
tween 10 and 70 GPa at room temperature [ |l4[|il| , a dis- 
crepancy that likely originates from the smaller O2-O2 
distance we calculate. Since the decline in intermolecu- 
lar distance associated with the Cmcm phase results in 
the formation of a permanent molecular polarization and 
intermolecular covalent bond, large dynamical effective 
charges (and hence absorptivity) would then be expected. 

The Raman- active vibron (A ff ) is in the range 1500- 
1650 cm , also in good agreement with measurement. 
Additionally, we predict five other Raman-active inter- 
molecular modes, and two can ostensibly be assigned 
to observed peaks below 500 cm -1 [jl5|,[53|, the B 3g and 
lowest-frequency B 2fl modes. Since our calculations over- 
estimate the IR-active mode frequencies, we regard this 
specific assignment as tentative. Extended scattering ge- 
ometries with fully polarized radiation may be needed 
to elucidate three further Raman-active modes predicted 
here whose intensities may be weak, given the normal 
mode displacement patterns associated with the linear 
structure and given the very anisotropic molecular po- 
larizability. In this context we note that the modes of 
interest are close in frequency to observed Raman peaks 
p3[ which have been assigned, so far, to possible over- 
tones and combinations of librons. It is also worth noting 
that even more modes would be expected of a more com- 
plicated e-phase, should it have a larger unit cell. Our 
results therefore suggest that an assessment of the phase 
diagram and lattice dynamics at low temperatures and 
over a wide pressure range will be of considerable exper- 
imental interest. 
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